Shortcut to adiabatic passage in two and three level atoms 
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We propose a method to transfer the population and control the state of two-level and three-level 
atoms speeding-up Adiabatic Passage techniques while keeping their robustness versus parameter 
variations. The method is based on supplementing the standard laser beam setup of Adiabatic 
Passage methods with auxiliary steering laser pulses of orthogonal polarization. This provides a 
shortcut to adiabaticity driving the system along the adiabatic path defined by the standard setup. 

PACS numbers: 32.80.Xx, 33.80.Be, 32.80.Qk, 03.65.Ge 



Introduction. — Two major routes for manipulating the 
state of a quantum system with interacting fields are 
based on resonant pulses or on adiabatic methods, such 
as "Rapid" Adiabatic Passage (RAP), Stimulated Raman 
Adiabatic Passage (STIRAP) and their many variants. 
In general terms simple fixed-area resonant pulses may 
be fast if intense enough, but quite unstable with respect 
to errors or fluctuations of the parameters, whereas adia- 
batic passage is robust but slow. For many applications, 
from Nuclear Magnetic Resonance (NMR) to quantum 
information processing, the ideal method should be fast 
and robust, combining the best of the two worlds. These 
two requirements are particularly demanding if quantum 
computing is to become feasible at all. It is possible 
to make the pulses more stable by combining them into 
pulse sequences, but in practice their use is limited by 
the longer times required with respect to the single pulse, 
the need to control phase angles and pulse durations ac- 
curately, or off-resonant excitations due to sharp pulse 
edges [1]. Moreover the error compensating properties 
of square-pulse sequences are not preserved when sub- 
stituting them with smooth pulses so that the design of 
good sequences requires "a good portion of experience 
and magic" 2]. In NMR, composite pulses are increas- 
ingly superseded by adiabatic passage methods [l| , which 
have also been very successful in chemical reaction dy- 
namics laser cooling Q, atom optics metrology 
@, interferometry Q, or cavity quantum electrodynam- 
ics [1, IU . When robustness is the primary concern, they 
are quite sufficient, and have as well become basic op- 
erations for quantum information processing, either to 
design rob ust g ates lfjl 11 | or in quantum adiabatic com- 
puting 
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13], which relies on an adiabatic evolution of 
the ground state from an initial to a final Hamiltonian. 
If speed is also important, however, the limitations may 
be severe [13j . Given the stated difficulties of compos- 
ite pulses, it is then quite natural to look for robustness 
and high operation velocities taking the adiabatic meth- 
ods as the starting point and shortening their duration 



somehow. Our objective here is to propose a shortcut to 
adiabatic passage (abbreviated as "SHAPE" hereafter) 
using a recent formulation by Berry |14| of "transition- 
less quantum driving" , related to work by Kato [HI on 
the adiabatic theorem. The specific applications we shall 
discuss are speeded-up versions of (2- level) RAP and (3- 
level) STIRAP schemes, as canonical examples of other 
adiabatic methods. Variants such as fractional RAP or 
STIRAP, and multilevel schemes may be treated along 
similar lines. 

The philosophy of the transitionless quantum driving 
algorithm [l3| is to supplement the Hamiltonian Ho(t) 
of a reference system with an auxiliary term Hi(t) to 
steer the dynamics exactly along the instantaneous eigen- 
states |A n (i)) of H {t) without transitions among them, 
formally in an arbitrarily short time, 



Hx(t) 



ihJ2(\ 9 t x n}{Xn\ - (A„|9 t A„)|A n )(A n |). (1) 



At variance with Lyapunov-control methods 16[ , the ex- 
tra term is independent of the time-dependent state so it 
leads to simpler, linear dynamics, and moreover it pro- 
vides systematically exact solutions for adiabatic follow- 
ing without the need for a trial and error approach to 
find a good control field [16J. Regarding its physical re- 
alizability, in general there is no guarantee that Hi may 
be easy to implement and each case needs a separate 
study. For example, for an Hq describing a particle in 
a time-dependent harmonic potential, Hi turns out to 
be a non-local interaction, and its realizability in a use- 
ful parameter domain for cold atoms remains an open 
question [13] ■ (Transitionless dynamics for the harmonic 
oscillator may be achieved with a local interaction by 
inverse engineering the frequency dependence with the 



aid of Ermakov-Riesenfeld invariants 18], or with state 
acceleration techniques t 19].) For a particle with spin 
in a time dependent magnetic field, Hi becomes a com- 
plementary, time-dependent magnetic field [3]. For the 
atomic two- and three-level systems studied here, Hi will 
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involve laser interactions added to the original laser setup 
implied by Hq as discussed below. 

Rapid Adiabatic Passage. — Let us consider first the 
speeding-up of a standard Rapid Adiabatic Passage that 
inverts the population of two-levels of an atom, |0) and 
|1), by sweeping the radiation through resonance. This 
broadspread technique originated in Nuclear Magnetic 
Resonance 20] but is used in virtually all fields where 
2-level systems may be controlled by external interac- 
tions, such as laser-chemistry, modern quantum optics or 
quantum information processing. When the frequency 
sweep is much shorter than the life-time for spontaneous 
emission and other relaxation times, it is termed rapid 
adiabatic passage (RAP). 

Using the rotating wave approximation, the Hamilto- 
nian Ho(t) in a laser- adapted interaction picture can be 
written as 

_ n ( A(t) n R (t) \ 
Ho{t >-2{n R (t) -A(t)J' (2) 

where £Ir(£) is the Rabi frequency, which we take to be 
real, and A(i) = usq—lul the detuning, assumed to change 
slowly on the scale of the optical period. It is the differ- 
ence between the Bohr transition frequency and the laser 
carrier frequency u>l, due to a change in the carrier fre- 
quency or a controlled alteration of the Bohr frequency 
by Zeeman or Stark shifts. The instantaneous eigenvec- 
tors are 

\X+(t)) = cos[0(t)/2]|l) - wx[0(t)/2]\O), (3) 
|A_(t)> = S in[0(t)/2]|l>+cos[0(i)/2]|O), (4) 

with the mixing angle 6(t) = arccos[— A(t)/fi(t)] 
and eigenvalues E±(t) = ±7ifi/2, where = 
yj A 2 {t) + £l R (i). If the adiabaticity condition 

i|fi Q | « |fi(t)|, (5) 

where fl a = [fl R (t)A(t) -fl R (t)A(t)}/fl 2 , is satisfied, the 
state evolving from \ip(t = 0)) = |A±(0)) follows the adi- 
abatic approximation 

|V±(t)>=exp|-ijf dt'E ± (t')}\\±(t)}, (6) 

whereas transitions will occur otherwise. Different adia- 
batic passage schemes correspond to different specifica- 
tions of n R and A for which ip± passes from one bare 
state to the other. The simplest one is the Landau-Zencr 
scheme with constant VL R and linear-in-time A. For the 
examples below we shall use the more adiabatic (and 
thus potentially faster) Allen- Eberly scheme [U l22j ]: 
tt R = ft sech(7rf/2^ ), A = {2/3 2 t /ir) tanh(7rf/2i ). Re- 
gardless of the scheme chosen Hi (t) takes the form 
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FIG. 1. (color online). (Dimensionless) energies e = E/(/3h) 
versus (dimensionless) time T = t/3 for the AE scheme: Di- 
abatic energies, dot-dashed (blue) lines; adiabatic energies, 
dashed (orange) lines; average energy evolving with Ho, solid 
(purple) line; average energy evolving with H — Ho + Hi, 
dotted (red) line, indistinguishable from the lower adiabatic 
energy, oj — 5, r = 1.22. 



where (up to a phase factor) fi a plays the role of the Rabi 
frequency of the auxiliary field. In principle H = Hq + Hi 
drives the dynamics along the 7?o-adiabatic path in ar- 
bitrarily short times, but there are practical limitations 
such as the laser power available. Moreover, a compari- 
son with i/o-dynamics is only fair if |O a | is smaller or ap- 
proximately equal to the peak Rabi frequency with the 
original laser setup. Independently of the scheme cho- 
sen and in a range of interaction times that break down 
the adiabaticity condition, it is remarkable that the dy- 
namics can be driven along the _ff -adiabatic path while 
fulfilling the inequalities |f2 a | < |f2| < |Oq|. 

The physical meaning and realizability of the auxiliary 
term are determined by going back to the Schrodinger 
picture: it represents a laser with the same time de- 
pendent frequency of the original one, but a differently 
shaped time-dependent intensity, and perpendicular po- 
larization. 

For the Allen-Eberly scheme the population of the ex- 
cited state Pi starting from the ground state depends on 
the dimensionless parameters r = to/3 and oj = Qq//3 [22j ]: 
Pi = l-sech 2 (2T 2 /7r)cos 2 [T(w 2 -4r 2 /7i- 2 ) 1/2 ]- A popula- 
tion transfer near to one (Pi > 0.999) and stable versus 
parameter variations is achieved for uj > 3 and t > 3. 
We may calculate tt a and the minimal time for which 
the maximum of Q a with respect to t is < fl - in the 
stated range this is accurately given by r m = tt/(4w), 
or to,m = 7r /(4^o)- The reduction factor with respect 
to the adiabatic time t q w 3 may be very significant, 
to,m/t a « tt/(12w), this is 0.09 for w = 3, or 0.01 for 
uj = 20. Of course the SHAPE Hamiltonian H may also 
drive the system along the adiabatic path outside the 
w, t > 3 domain as illustrated in Fig. [T] 

For comparison, the population of the excited state due 
to a square 7r pulse with on-resonance Rabi frequency flo 
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(8) 



Complete population transfer requires A = 0, and a pulse 
time ffl = jf- • For the same Qq and limiting the auxiliary 
laser by |f2 a f < fig, the minimal characteristic time to,m 
of the SHAPE method is of the order of ta, io,m = *r/4. 
In fact the actual interaction time to implement a suc- 
cessful population inversion with the AE scheme (SHAPE 
corrected or not) should be a few times to] this may be 
estimated from the dependence of the excited popula- 
tion of the adiabatic state with time, which is > .999 for 
t > 8t . 

Figures 2 and 3 show examples of the fidelity (Pi) 
with respect to variations in the Rabi frequency and 
detuning with SHAPE (AE scheme for Ho), the evolu- 
tion with H (AE scheme), a Rabi 7r-pulse, and a com- 
posite ^(x)7r(y)^(x) pulse, a fault-tolerant combination 
where x, y refer to the laser polarization (and Pauli ma- 
trix) involved. Clearly SHAPE provides a fast, robust 
and efficient population inversion compared to all other 
methods. All cases are for the same fio, and in SHAPE 
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FIG. 2. (color online). Dependence of the fidelity (Pi) on 
the changes of (a) the Rabi frequency from fio to fio(l + n) 
(b): the detuning by S. SHAPE, AE scheme (solid, red 
line) with il a < fio; ordinary Adiabatic Passage (dotted, 
black line); Rabi pulse (dashed, blue line); Composite pulse 
~(x)w(y)Z(x) (dot-dashed, purple line), Qo = 2ir x 5 MHz; 
t = 25 ns; p = 2tt MHz. 

Stimulated Raman Adiabatic Passage. — Similar ideas 
can be applied to three-level STIRAP. The Hamiltonian 
H (t) for the two-photons resonance case within the ro- 
tating wave approximation (RWA) and in a laser adapted 
interaction picture reads Q 

( o n p (t) o 

H (t) = - n p (t) 2 A fi a (t) 

1 \ o n s (t) o 

in terms of the Rabi frequencies for the Stokes, f2 s (t), 
and pumping lasers, O p (t), and the laser detuning hA = 
(E 2 — Pi) — hojp = (E 2 — E 3 ) — huj s . The instantaneous 
eigenstates |A n ) are 

\X+(t)) = sin6»sin0|l) + cos0|2) + cos 9 sin 0|3), 

| A_ (t)) = sin(9cos0|l) — sin<^|2) + cos cos 0|3), 

|A (*)> =cos6»|l) — sin6>|3), (10) 



(9) 



with eigenvalues given by E + (t) = UClcot((j>/2), E = 0, 
and E-(t) = — 7ifitan(</>/2). The time-dependent mix- 
ing angles 9 and 4> are respectively defined by tan# = 
n p (t)/Q s (t) and tan(20) = Q/A(t), whereas fl = 

The population transfer 1 -> 3 is re- 



(11) 



n|(t) + n2(t). 

alized by the "dark state" |Ao). 

The Hamiltonian Hi(t), takes now the form 



Hi(t) = in 



4> sin t 
\ sin 9 
— 6 cost 



'cos 9 




with 9 = [ti p (t)tt s (t) - n s (t)fl p (t)]/n 2 , and 

= [(ti p {t)n p (t) + n s (t)n s (t))A(t)]/[2n(A 2 + n 2 )]. 

We would need in principle three new lasers to im- 
plement this Hamiltonian. The ones connecting levels 
1-2 and 2-3 should have the same frequency as the 
original ones but orthogonal polarization, and the field 
connecting levels 1-3 should be on resonance with this 
transition to get an interaction picture Hamiltonian like 
(TTTI) (The RWA approximation is assumed in all cases.) 
If 1-3 is electric-dipole-forbidden, a magnetic dipole 
transition may be used instead. If we are only interested 
in performing a full passage from 1 to 3 and do not 
want to reproduce all the effects of the full Hamiltonian 
Hq + Hi, Hi may be simplified by retaining just the 1-3 
interaction, 



H[(t) = 




(12) 



where fl' a = 29. That this is so may be seen by work- 
ing out the Schrddinger equation in the adiabatic basis: 
d(Xo(t)\tp(t))/dt does not depend on <j) so that the 1-3 
and 2-3 auxiliary lasers can be left out without affecting 
(Xo(t)\ip(t)) . In the examples below we have chosen the 
pulse shapes Q 



n p (t) = Cl (t)f(t-T); 

fit) - { f {nt/T) 



n,(t) = n (t)/(t), 

(0 < t < T), 
(otherwise). 



(13) 



Fig. [3] shows a STIRAP Stokes-Pump pulse sequence 
where T is too short for complete population transfer 
because adiabaticity breaks down, see Fig. @h.. We can 
remedy that with the auxiliary interaction in (|12[) . see 
Figs. [3] and 2b. Keeping \Cl' a \ < £l the process duration 
is reduced approximately ten times with respect to the 
ordinary STIRAP scheme. 

Discussion and conclusions. — A method to achieve fast 
and robust population transfer in two-level and three- 
level atomic systems has been presented, based on sup- 
plementing the laser setup of standard adiabatic passage 
methods (RAP or STIRAP) by additional, properly time- 
shaped pulses with orthogonal polarization. The Hamil- 
tonian Pi (t) that describes the additional steering pulses 
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FIG. 3. (color online). Time evolution of Rabi frequencies in 
a STIRAP sequence of laser pulses defined by Eq. (|13[) with 
n = 2tt x 5 MHz, A = 2tt x 0.5 MHz, and T = 0.26 /js. 
Dashed blue line: fl s /fio; Solid red line: Q p /Qo; The dotted 
orange curve represents Q' a /Qo. 
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FIG. 4. (color online). Time evolution of the populations 
of levels 1, 2 and 3 for STIRAP: (a) Hamiltonian H {t); (b) 
Hamiltonian H(t) with additional H\(t). Same parameters as 
in Fig. El 



providing a shortcut to adiabaticity is given by a general 
algorithm to drive quantum systems without transitions 
[14j . Other states (such as superpositions of two-levels) 
may be prepared by speeded-up (SHAPE) versions of 
fractional RAP or STIRAP, and, if necessary, the phases 
may be controlled thanks to the freedom to choose the 
reference Hamiltonian Ho(t) or phase gates. As an out- 
look, similar techniques may provide a way to carry out 
adiabatic computation in a finite time 3, 25 1, or to speed 
up logic gates based on adiabatic processes [2|| 27 [, in- 
terferometric techniques in superconducting qubits [28[ 
or quantum dots [29[ , and the creation of entangled pairs 
of two-state systems [1(3]. The SHAPE method is com- 
patible with approaches that optimize Hq such as the 
quantum brachistochrone approach [311 ] . since, after op- 
timizing the adiabatic process, it leads to the design of 
even faster process. Other adiabatic techniques [32[ may 
also benefit from these ideas speeding up the avoided 
level crossings and keeping their stability versus param- 
eter variations. 
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